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We investigate the validity of the approximate method to describe a strong gravita- 
tional lensing which was extended by Alard on the basis of a perturbative approach to 
an Einstein ring. Adopting an elliptical Navarro-Frenk- White (NFW) lens model, we 
demonstrate how the approximate method works, focusing on the shape of the image, 
the magnification, caustics, and the critical line. Simplicity of the approximate method 
enables us to investigate the lensing phenomena in an analytic way. We derive simple 
approximate formulas which characterise a lens system near the Einstein ring. 
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1. Introduction 

Cold dark matter is one of the most important components in the universe. The cos- 
mic microwave background anisotropics and the large scale distribution of galaxies 
cannot be naturally explained without the cold dark matter component. The mean 
density parameter of the cold dark matter has been measured precisely,^ but its true 
character has not been identified. The elementary particle physics predicts possible 
candidates of the cold dark matter, and many experiments are ongoing aiming at a 
direct measurement. 

The cold dark matter is considered to be distributed associated with each galaxy, 
forming dark matter halo. Then, the investigation of the structure of the halos is 
quite important in exploring the nature and the origin of the dark matter. The 
strong gravitational lensing is a useful probe of the halo-structure (see, e.g. ,121 for a 
review). Especially, a lens system near Einstein ring is useful because a wealth of 
information can be obtained.^ 

Besides, the strong lensing systems are also useful as a tool of the dark energy 
study.SHSl Because of the recent observational developments, many strong lensing 
systems have been found. The strong lensing statistics is now becoming one of the 
powerful tool for exploring the nature of the dark energy.!^ Future dark energy 
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surveys will detect much more strong lensing systems (see, e.g.,tl^, and the strong 
lensing system will play a more important roll in cosmology. 

In realistic situations, the mass distribution in a halo is not simple, which makes 
reconstruction of the lens model complicated. The lens equation is complicated for 
a non-spherical lens model, which needs to be solved numerically. Then, analytic 
approximate approach to strong lensing system is useful, if its validity and accu- 
racy are guaranteed. A perturbative approac h to th e lensing system close to the 
Einstein ring configuration was developed, e.g.j ^^^^^ . Recently, Alard extended the 
perturbative approach, which is applied to analyse lensing systems IISHISI 

In the present paper, we investigate the validity of the perturbative approach to 
the lensing system close to an Einstein ring, assuming an elliptical lens model. We 
demonstrate the validity of the perturbative approach quantitatively, by comparing 
with an exact approach on the basis of the numerical method, focusing on the 
shape of the image, the magnification, the caustics, and the critical line. Using 
the approximate method, expanded in terms of the ellipticity parameter of the 
lens model, we derive simple approximate formulas which characterise an elliptical 
lensing system near the Einstein ring in an analytic way. 

This paper is organized as follows: In section 2, we briefly review the basic for- 
mulas for the gravitational lensing and the per turbative approach to a perturbed 
Einstein ring, based on the work by Alard.ll^In section 3, we compare the pertur- 
bative approach with the exact approach that relies on a numerical method, focusing 
on the shape of lensed images, the caustics, the critical curve, and the magnifica- 
tion, respectively. We demonstrate the validity of the perturbative approach at a 
quantitative level. In section 4, some useful formulas are presented, which are de- 
rived using the perturbative approach in the analytic manner. Section 5 is devoted 
to summary and conclusions. Throughout the paper, we use the unit in which the 
speed of light equals 1. 

2. Basic Formulas 
2.1. General basis 

We briefly review basic formulas for the strong lensing (e.g.^^. The deflection angle 
of a lens object is determined by 



where $ is the gravitational potential of the lens object, x is the radial coordinate 
connecting the observer and the lens object, ^ is the two dimensional vector on the 
lens plane, which is orthogonal to the coordinate x- The gravitational potential $ is 



where A denotes the 3-dimensional Laplacian, and G is the gravitational constant. 




(1) 
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Introducing the surface mass density S(^), which is the projected mass density 
on the lens plane, 

m= P{i,x)dx, (3) 
Jo 

and the lensing potential, 

m= H€x)dx, (4) 
Jo 

which are related by 

A(2)V;(a = 47rGS(^"), (5) 
where A^^' denotes the 2-dimensional Laplacian. The solution is 

^(1) = 2G'y (f^'T,{i') log ll*- I*'! + constant. (6) 
Then, the deflection angle is 

& = 2^$) =4G f dH'm')4^. (7) 
Now we consider the gravitational lens equation, 

V=^i-Dj.s&, (8) 

where Dls is the angular diameter distance between the lens and a source object, 
Ds is the distance between the observer and the source, I?l is the distance between 
the observer and the lens, and ff is the two dimensional vector on the source plane 
orthogonal to the coordinate x- Introducing a characteristic length in the lens plane, 
^0, and ?7o = ^qDs/Dl in the source plane, we define 

y = v/vo (9) 

x = i/^o, (10) 

then, the lens equation becomes 

y = x-\/s(l>{x), (11) 



where we defined 



with 



(f) = (fx'K{x')log\x-x'\ (12) 



K{x) = (13) 

^cr 

S^=^^«. (14) 
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2.2. Perturbative approach to Einstein ring 

Next, we review the perturbative approach to the Einstein ring developed by Alard 

115116J unni) 

When the projected density of the lens S is circularly symmetric 
and the source is located at the origin of the source plane, y = 0, an Einstein ring 
is formed. The radius of the Einstein ring is determined by 

where (/)o(|af|) denotes the circularly symmetric lens potential. We denote the solu- 
tion of Eq. (1151) hy X — Xe- Thus, \xe\ is the Einstein radius. Hereafter, we use the 
notation x = \x\ and xe — \xe\- 

We consider the perturbative approach to the Einstein ring, then assume that 
the deviation from the circularly symmetric lens is small. Introducing the small 
deviation, which are denoted by the quantities with S, 

y^Sy, (16) 

X — Xe + Sx, (17) 

(j){x) = Mx)+SH^), (18) 

the lens equation ([TT]) is rephrased as 

Sy = XE + 5x -\/s[(l)o{x) + 6(l){x)] ^ ^ ^. (19) 

x—xe+Sx 

Assuming that the deviation from the circularly symmetric lens is small, we intro- 
duce the small expansion parameter e. Explicitly, we assume 

Sy^O{e), (20) 

Sx = 0{e), (21) 

(50(f) 0{e). (22) 
We find the following lens equation at the lowest order of e, 

Sy^6x-[{6x-\/s)yxMx)+'^xS(j){x)] ^ ^ , (23) 



where we used Eq. (flSl). 

We consider a circular source with the radius Svs , whose centre is located at the 
coordinate {Syio,Sy2o) on the source plane. Then, the circumferences of the source 
is parameterised as 

Syi = Syio + Srs cos <f (24) 
Sy2 ^ Sy2o + Srs sin (fi (25) 

with the parameter Lp in the range < ip < 27r, where we assume that 6yio, Sy2o 
and Srs are the quantities of the order of e. Similarly, we may rewrite the image 
position X = Xe + 5x sls 

xi = {xe + Sx) cosO, (26) 
X2 — {xe + Sx) sin6. (27) 
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Here, we don't need to co nsid er the perturbation of 9 because of the symmetry 
of the un-perturbed imageff^ In the un-perturbed situation, the image of a point 
source is a circle, and there are an infinite number of the image at any 9. Then, at 
any angular position of the perturbed point, there is always an un-perturbed point 
at the same angular position on the circle. 
The lens equation p3|) yields 

Syio + SrsCosf = [Sx cos9 — dl(j)o{x)) 
sinC 



cos9dx64>{x, 9) 



-de5(t>{x,9)\ 

Ju 

^Vin ^ sSvciip — [(5a;sin6' (l — 9^(/)o(a:)) 
— sin6'9j;(50(a;, 6') — ^^^^ 



-de8<\>{x,9)\ 



Combining these equations, we have 
1 



bx = 



where we defined 



yed(p(x, t 



Syio sin 9 + 6y2o cos 9 



(28) 
(29) 



dxS(p{x, 9) + 5yio cos 9 + (5y2o sin 9±^/^?{x,9) 



(30) 
(31) 



This is the formula derived by AlardP^GH 

For comparison, we summarised the corresponding formula without the pertur- 
bative approximation in the appendix A. One can derive Eq. (j30p from Eq. (jA.Tp . 



2.3. Perturbative NFW lens model 

A mass model commonly used for strong lensing is based on high-resolution numer- 
ical simulations of dark-matter halos in the A CDM framework by Navarro, Frenk 
and White (1996, 1997; hereafter NFV^EG^l)^ 

in which the density profile is pa- 
rameterised by the scale radius and the constant p^, 

Ps 



p{R) 



{R/rs){l + R/rs) 



2 ' 



(32) 



where R = Vs+X^ is the 3-dimensional length. Choosing Ts = fo, the lens po- 
tential becomes 



0NFw(a;) = ^i^sF{x), 



where we defined 



F(a;) = iw| + <( 



h2arctan2 ^ ^ [x > \) 

Y x-\-i ^ ' 

~2arctanh^./T^ (x <\) 
logx {x = 1) 



(33) 



(34) 
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and 

= (35) 

We denote the solution of the lens equation for the circular NFW lens model by uq, 
which satisfies 

d 

Uq = 4>^Yw{uq). (36) 

In the present paper, for an asymmetric lens model, we adopt the potential 

(l){x) — 0NFW {x\Jl — ?;cos 26^ . (37) 

Instead of x and y, we introduce x and which is normalised by the Einstein radius 
Wo defined with Eq. ((36|) . 

X = uqx, (38) 
y = uoy. (39) 

Then, the lens equation is rewritten 

^ = J - V.^ (40) 

with 

^=^F (uoxJl~f]cos 29] , (41) 
where x = |a;|, F{x) is defined by Eq. and 



= log y + 



'"0 


-1 


"0 


+1 


/I- 


Uq 



(42) 



0^ V 

1 (Uo = 1) 

Finally, the potential of the elliptical NFW lens is written as 

4>^M^) + T4>{x,e) (43) 

with 

M^) = ^F{uox), (44) 

S^{x, 0) = ^ (F{uoI^/l-r|cos2e) ~ Fiu^x)] . (45) 

In the appendixes A and B, useful formulas related with the elliptical NFW lens 
potential are summarized. In these appendixes, only the case uqX\J\ — 77 cos 29 < 1 
is described, but the case uox^/1 — 77 cos 29 > 1 is obtained by the analytic contin- 
uation. 
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3. Validity of the Perturbative approach 

We here investigate the validity of the perturbative approach, comparing with re- 
sults without any approximation. For being definite, we consider the following three 
cases. 

(a) Exact approach without any approximation. 

(b) Perturbative approach described in the previous section. 

(c) Approximate approach: the perturbative approach (b) plus the lowest-order 
expansion of 77 (See also the appendix C). 



3.1. Image 

We consider the lensed image of the circumference of the circular source, whose 
center is located at {Syio, Sy2o)- The source's radius is Sr^. In the exact approach 
(a), the circumference of the lensed image is obtained by solving Eq. (jA.7p . which 
is derived in the appendix A. Eq. (jA.7|) can be solved with an iterative method 
numerically. On the other hand, in the perturbative approach (b), the circumference 
follows Eq. which is equivalent to 



6x = 



+ Syio cos 9 + dy2o sin 6 ± w (cc 



(46) 



with 



A^{x,9) = drj^ - (^-^d0S(j)(x,e) - 5yiosme + Sy2aCOS0^ . (47) 

Eqs. dH]) and (gT]) are the same as Eq. (12) in the reference Alard In the 
approximate approach (b), we use Eqs. (jB.l| ~ (jB.3P to find the solution, Eq. (|46|) . 
In the approximate approach (c), we use the quantities of the lowest-order of 77, 
given by Eqs. (fC^ and ([03]) . instead of Eqs. ([R2l) and (|R3| . 

The upper left panels of Fig. [1] show the lensed image, (a) is the exact approach, 
(b) is the perturbative approach, and (c) is the approximate approach, respectively. 
The panel (d) plots these three approaches for comparison. Here we adopted the 
parameters 5rs = 0.07, SyiQ = 0.09, 5y2o — 0, rj — 0.15 and uq = 0.5. In this case, 
the lens effect splits the image into four. The dashed circle in each panel is the 
Einstein radius of a point source. 



It is useful to s umm arize the differences of the notation between the present paper and the 
reference by Alard CSl. Our equa tion s can be obtained by the following transformation from the 
equations in_jthe paper by Alard El, r — >• ir, dr ^ 5x, 4>0 — >• <l>0, V' ~^ ^4'i Ro ~^ ^I's, — >■ 5j/io, 
and j/o <5j/20- 
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Fig. 1. Comparison of the three approaches, (a) exact approach, (b) pcrturbativc approach and 
(c) approximate approach of the lowest-order of the expansion of rj. The panel labelled by (a), (b) 
and (c) corresponds to the approach (a), (b) and (c), respectively, while (d) plots all the three 
approaches for comparison. The upper left panels show the lensed image of the circumference of 
a circular source. The dotted circle is the Einstein ring. The long and short dashed curve is the 
critical line. The parameters Srs = 0.07, Syio = 0.09, <5y20 = Oi ^ = 0.15 and uq = 0.5. This is 
type I. The upper right panels show the magnification factor fi on the circumference of the images 
as a function of the angular coordinate of the image plane. The lower left panels show the critical 
line. The dotted circles are the Einstein ring. The lower right panels show the caustics on the lens 
plane. The circle is the circumference of source. 
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3.2. Magnification factor 

The magnification factor due to the gravitational lensing is given by the inverse of 



the determinant of the Jacobian matrix. 



fi=\J\- 



12 



ox 



(48) 



For a general lens potential, the determinant of the Jacobian matrix can be written 
as 



X 



V dx^ / V dx X 89^ / x\x 89 



18^ 



1 /I 



8^ 



8x89J ' (^9) 

where we used Eqs. (jA.10[) ^ (jA.14l) . This is the exact approach (a) for the magnifi- 
cation. 

On the other hand, in the case of the perturbative approach (b), we have 



J 



1 



& 

8x'^ 



dx — 5x 



92 



92 



dx^ dx 



89^ 



(50) 



where we used the condition of the Einstein ring. 8(^qI8x\x=\ = 1. Note that 
0(x,6') = 4>o{x) + 5(j){x,9) and x = 1 + 6x. The formulas (|B.1|) ^ (|B.4|) are used 
for the perturbative approach (b), and (jC.2[) ^ (|C.4[) for the approximate approach 
(c), respectively. 

The upper right panels (a)^(d) of Fig. [Ijshow the magnification of the circum- 
ference of lensed image as a function of the angular coordinate of the image plane. 
The panels (a)~(c) correspond to the three approaches (a)~(c), respectively. The 
panel (d) plots the three approaches. The model parameters are the same as those 
of the upper left panels for the lensed image. 



3.3. Critical line 

The definition of the critical curve is J = on the image plane. In the exact 
approach (a), we solve 

la^ 1/ &^Y^/ld^ 
di^^W^^y ~"ds!^) [^'89 ~ 8i89J ' ^ ' 

which is obtained using Eq. (j49|) . 

In the perturbative approach (b), the critical line is 

~ _ 1 (8T4> 8^H 

^ ^ 1 - 8^^0/di^ \ 

from J = with Eq. (|50p . This equation is the same as Eq. (30) in the reference by 
Alard^^. In the perturbative approach (b), we use Eqs. (|B.2|) and (jB.4p . while we 
use Eqs. (jC.2p and (jC.4[) in the approximate approach (c). 

The lower left panels (a)^(d) of Fig. [T] show the critical line (solid curve) on 
the image plane. The panels (a)^(c) correspond to the three approaches (a)~(c). 
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respectively, and the panel (d) plots the three approaches. The dotted circle is the 
Einstein ring. The model parameters are the same as those of the other panels of 
Fig.ffl 



3.4. Caustics 

The caustics are defined by J = on the source plane, which can be mapped from 
the critical line on the image plane by the lens equation. The caustics are important 
in understanding the nature of deformation of an Einstein ring. 

In the exact approach (a), the caustics are obtained by substituting the solution 
of Eq. dsn into (|X5|) and fO]) . 

In the perturbative approach (b), the caustics are given by 



5V2 



09^ 



x = l 



sin 9 



x = l 



89 
~d9 



sin6', 



(53) 
(54) 



x = l 



which are obtained by substituting Eg. ([5^ into the lens equation ([^^ and 
Eqs. (|53|) and (f54| are the same as Eq. (31) in the reference by Alarcjl^, where Xs 
and ys are used instead of 5yi and 5y2, respectively. In the approximate approach 
(c), we use ([C3t and (Ia4| . 

The lower right panels (a)^(d) of Fig. [T] show the caustics (solid curves) on 
the source plane. The panels (a)~(c) correspond to the three approaches (a)^(c), 
respectively, and the panel (d) plots all three approaches. The spherical circles are 
the circumference of the source. The model parameters are the same as those of the 
other panels of Fig. [T] 



3.5. Typical configurations 

In Fig. [H we adopted the parameters 5rs = 0.07, 5yiQ = 0.09, 5y2o = 0, V — 0-15 
and Mo = 0.5. Figs. [2] and [3] are the same as Fig. [1] but with rj — 0.1274 and 
T] = 0.08, respectively, instead of ?7 = 0.15. The other parameters are the same as 
those of Fig. [1] then the source configuration of Figs. [1] [21 and |3] are the same. 
As T] becomes smaller, the size of the caustics becomes smaller, and the shape of 
the critical line becomes more spherical. Also, as rj becomes smaller, the right-side 
three separated images become to merge. Fig. [2] is the critical configuration when 
the merger occurs. One can observe that the merger occurs when the circumference 
of the source contacts with the caustics. In the upper right panels of Figs. [2] and [U 
the large enhancement of the magnification appears. This refiects the facts that the 
magnification diverges when a source is on the caustics and when the image crosses 
the critical line. Figs. [TJ [2] and [3] represent typical types of lensed image, which we 
call type I, II and III, respectively. 
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As T] becomes smaller furthermore, the left-side image and the right-side image 
become elongated, and form a ring, as is demonstrated in Figs. |4] and [5l The pa- 
rameters of Fig. |4] are 5rs = 0.07, 5yia = 0.06, 5y2Q — 0, — 0.0205 and uq = 0.5. 
Fig. [5] is the same as Fig. |4]but with the different value of — 0.011, instead of 
r] = 0.0205. These two images are almost rings, i.e., Einstein rings with finite width. 
Fig. [4] is the critical configuration when a ring is formed. Figs. |4] and [5] represent 
typical types of lensed image of a ring, which we call type IV and V, respectively. In 
these cases, the size of the caustics is smaller than the source size, which is clearly 
shown in Figs. 0] and [SJ We note that the critical configuration, type IV, appears 
when the circumference of the source comes in contact with the caustics at the left- 
side caustic. The divergence of the magnification appears when the source overlaps 
the caustics and when the image crosses the critical line. 

3.6. Validity of the perturbative approximation 

Let us discuss about the validity of the perturbative approach comparing with the 
exact approach. Fig.lHlplots the outer position of the image at 6* = 0, i.e., x{6 = 0), as 
a function of 77. (See also Fig. [71) The curves labelled by (a), (b) and (c) correspond 
to the three approaches. In the exact approach (a), x{9 = 0) is given by solving 
Eq. (IA.7I) . In the perturbative approach (b) and the approximate approach (c), 
x{9 = 0) = 1 + Sx{9 = 0), where Sx{6 — 0) is given from Eq. pS)) with the 
formulas in the appendix B and in the appendix C, respectively. We are considering 
the circumference of a circular source, and x{9 = 0) has the two solutions, which 
correspond to the two solutions of the sign ± in Eq. (jA.7[) in the exact approach 
(a) and Eq. (|46)l in the approaches (b) and (c). Fig. [6] plots the outer solution, 
which corresponds to the solution of the sign -I- in Eqs. (jA.7[) and p6)) . In the limit, 
77 — 0, (b) is equivalent to (c). The difference between (a) and (b) (or (c)) in the 
limit 77 0, comes from a limitation of the perturbative method. As one can see 
from Fig. [6l the difference between (a) and (b) becomes large for rj ^ 0.5, where the 
perturbative approach bre aks down. The similar result was obtained in Figure 1 in 
reference by Alard (2007) Besides comparing x{9 = 0), it is useful to compare 
the position of the image edge as examined in Fig. 8 in the below (See also the 
reference by Peirani, et al.(2008)^^). 

Next, let us focus on the critical configurations where the number of the images 
changes. Four images of the type 1 change to two images of the type III as rj becomes 
smaller. The type II is the critical configuration. The two images of the type III 
changes to a ring configuration of the type V as 77 becomes smaller. The type IV 
is also the critical configuration. Fig. [8] examines these critical behaviours. Fig. [8] 
shows the angles of the edge of the image, 9i, 6*2 and 9^, which are defined using 
Fig. [71 as a function of rj. In each panel, we plot 0i, ^2 and ^3, and the top curve 
is ^3. As 77 becomes small, 92 and 9i merge, when the critical configuration (type 
II) appears. We refer rju as the critical value of 77 when 02 and 9i merge. For 
77 < rju, we have no solution for 9i and ^2, In the upper left panel of Fig. [8l the 
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vertical dashed line labelled by I, II and III, are the value of 7], adopted in Fig. [1] [2] 
and [3l respectively. The solid curve is the exact approach (a), the dashed curve is 
the perturbative approach (b), and the long dash-dotted curve is the approximate 
approach (c). Note that the three approaches (a), (b) and (c) agree very well. 

The upper right panel of Fig. [8] is the same as the upper left panel, but with 
Syio = 0.06, instead of Syio = 0.09. Similarly, the lower left (right) panel is the same 
as the upper left panel, but with Syio — 0.15 ((5j/io — 0.3). Thus, the agreement 
between the three approaches is better for smaller values of Sy iq. Note also that the 
critical value yyn becomes larger as Syio becomes larger. 

In the upper right panel of Fig. [51 the vertical dashed line labelled by IV and V, 
is the value of 77, adopted in Figs. |3]and [SJ respectively. In this panel, 9^ increases 




Fig. 2. Same as Fig. [T]but with the parameters Sra = 0.07, Syio = 0.09, <5j/20 = 0, = 0.1274 
and Mo = 0.5. This is type II. 
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-1 1 -1 1 -0.3 0.3 -0.3 0.3 



Fig. 3. Same as Fig. [T] but with the parameters Srs = 0.07, Syio = 0.09, Sy20 = 0, rj = 0.08 and 
tto = 0.5. This is type III. 



as T] becomes small. We define rjiy by the smallest value of r/ with which we can 
find the solution of ^3. The type IV critical configuration appears at ?7 = 771V- as 
shown by the vertical dashed line labelled by IV. The critical configuration of type 
IV appears only in the upper right panel of Fig. |8l where a ring image appears for 
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4. Discussions 

As is demonstrated in the above, the approximate approach (c) is not so bad. 
An advantage of the approximate approach (c) is that the simphcity enables us 
to investigate the lensing phenomena in an analytic way, as we will show in this 
section. 



4.1. Condition of the critical configurations 

First, we consider the condition that the critical configuration of type II appears, 
which is the transition point that the number of the images changes from four to 
two. 
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In the perturbative approach, the condition {xe,9) > must be satisfied for 
the existence of the solution of the lens equation (see Eqs. (150]) and (PTjl ). In the 
approximate approach (c), from Eqs. (H71) and (|C.3p . we have 

A^{x = 1,6) = (5^/ - (?7 sin 26* - (^yTo sin 61)2, (55) 

where we used Sy2a = 0- Fig. IH] plots A'^{x = 1,6') as a function of 6 for the five 
typical cases, I ^ V, corresponding to Figs. [1] ^ [5j Note that our models satisfy 
A2(J = 1, 6*) = A2(5 = 1, 27r - 6*). Fig. H only plots the range of < 6* < tt. Lensed 
image appears under the condition A'^{x = 1,9) > 0. The number of separated 
regions with A'^{x = 1,6) > determines the number of the images. Also the zero 
points oi A^{x = 1,9) = m the panel (1) of Fig. [9] correspond to 9i, 6*2 and 6*3 in 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

7] 7} 

Fig. 6. Comparison of the solution of the lens equation x{6 = 0) as a function of rj. The solid 
curve is the exact approach (a), the dashed curve is the perturbative approach (b), and the long 
dash-dotted curve is the approximate approach (c). The upper left panel adopted Svs = 0.07, 
^yw = 0.09, 5y2o = and no = 0.5. The upper right panel is the same as the upper left panel but 
with the different value of 5yio = 0.09. The lower left (right) is Sts = 0.2 5j/io = 0.3, no = 0-3 
(no = 0.9). 



Fig.m 

Let us first explain the local minimum and maximum points of A^(a; = 1,9), 
9niin± and ^max, respectively, (see the upper left panel of Fig. [9]). From (|55|) . we 
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x(e=o) 



have, 




Fig. 7. Definition of the angles 6i, 62, and 63, for the image of type I. 



— A2(i =1,9) = sin 6l[-16?7^cos3 61 + 12?7(5yio 008^6' 
dO 

+ {STf - 25yio ) cos 61 - A-qSyw] 



= — 16?7 sinO cos( 



5yi 



277 



5yw + V^V_+^yio_ 
8?7 



X COS! 



5yw - V i2rf + Sy^ 



(56) 



We find that A^(a; = 1,0) has a local maximum in the range of < ^^ < tt at 



cos 9„ 



Syio 
2ri ' 



(57) 
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mill viv 




Fig. 8. 6i, 02 and 63 as a function of rj. In each panel, the top curve is 63. As r) becomes small, 
9\ and 62 merge, when the critical configuration (type II) appears. The solid curve is the exact 
approach (a), the dashed curve is the perturbative approach (b), and the long dash-dotted curve is 
the approximate approach (c). The dashed vertical lines are the values of -q adopted in Figs.[T)^{5l as 
labelled in this figure. In this figure v^e adopt Sts = 0.07, Sy20 = and no = 0.5, and Syio = 0.09 
(upper left panel); 5y\o = 0.06 (upper right panel); (5yio = 0.15 (lower left panel); Syio = 0.3 
(lower right panel). 



and the two local minimum at 



cos ^niini = ^ (syw ± ^^Iff+^y^ ] ■ (58) 



877 

The values of the local maximum and minimum are 



A2(J= 1,0 = 0,„,,) (59) 
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Fig. 9. A'^(x = 1,9) as a function of 9 (see Eq. lj55|l V Here the parameters of each panel (I) 
(V) are the same as those of Figs.[T]~ (5] respectively. 



and 



1 



^inini ) 



128772 



± 32riHyioJ 32if + 5yw + 128<5r,%2 



A 



±1 



(60) 



respectively. 

Furthermore, A2(-e — 1,9) also has other local maximum at 9 — Q and 9 = tt. 
From ((55l) . we have 



A2(x = 1,6* = 0) = A2(S = 1,6* = tt) = (5r 



(61) 



This means that the thickness of the arcs is approximately determined by Sr^ , which 
can be seen in Figs. [T}^ [5l 

Note that the critical configuration type II (IV) appears when A^ = (A^^ = 0). 
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From Eq. (|60p . the condition Aj_ = is rephrased as 

2 



/±(^)=_8(^j , (62) 



where we defined 



and ^ = 5yia/2ri. Fig. [10] plots /"^(^) and as a function of ^. Note that 

/±(0) = -8. From Fig. [TOl one finds > /"(O- This means that A^. > A?., 

which can be proved explicitly. This also means that the configuration type always 
changes as V — > IV — > III — > II — > I, as Syio/rj changes from infinity to 0. 




_60 I ' \ ' \ ' \ ' \ ' 1 

0.2 0.4 0.6 O.B 1 

Fig. 10. f'^ (upper curve) and /~ (lower curve) as a function of 

We can easily find the critical condition that the type II appears in an analytic 
way, as follows. The condition is A^ = 0, i.e., Eq. (|62|) of the + sign, which has a 
solution around ^ ~ 1 for rj Sy iq. By expanding /^(C) around ^ = 1, we have 

/+(c)^^(c-i)3 + o(e-ir- (64) 

With this approximation, A^ = yields 




2, ^ • 

Fig. [TT] is the diagram to show where the typical configurations appear on the 
r] — 5yiQ plane. On the boundary between I and III, the critical configuration type II 
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appears, while the type IV appears on the boundary between III and V. The dashed 
curves satisfy = 0, the dot-dashed curves satisfy A?_ = 0, which are obtained 
by solving Eq. ((62|) of the sign + and — , respectively. The solid curves plot Eq. (|65l) . 
The agreement between the dashed curve and the solid curve means the validity of 
the approximate formula of (|65|) . Here, we adopted Sr^ = 0.07, uq = 0.5 (upper-left 
panel) and 5rs = 0.2, uq = 0.3 (upper-right panel), respectively. The lower left 
(right) panel assumes mq — 0.5 (mq = 0.9). 

The critical boundary of Fig. [TT] is obtained on the basis of the approximate 
approach of the lowest order expansion in terms of -q. The exact critical value should 
be found by using figures like Fig. |8l The points in Fig. [11] are obtained by making 
figure like Fig. [51 The cross, the triangle, and the square are the results with the 
approach (a), (b) and (c), respectively. Thus, the critical boundary is estimated 
lower when we use the approximate approach, as shown in each panel of Fig. ITT] 
However, this figure also demonstrates that the approximate approach is quite good 
as long as ?7 ^ 0.3. 

The condition that the critical configuration II appears is = 0. We may 
write the condition as 

5yw + ri < 5rs. (66) 



4.2. Relation with caustics and critical line 



We consider the condition that the circumference of the source comes in contact 
with the caustics. The intersection point is obtained by substituting Eqs. ()53|) and 
(|54| into Eqs. ([24]) and ([25]), which yields 

2 / \ 2 

X—1 , 



(9612 



de 



Syio 



2Sy 



92 



10 



COS! 



dS<j) 
~d6 



sin0 = 5r 



(67) 



where we used 5y2o = 0. Within the approximate approach (c), we use Eqs. (|C.3p 
and (|C.4p . which we substitute into Eq. ([ST]) , then 

g{e) = 12772 cos'' e - 4Tj6jno cos^ - Uif cos^ 61 + 47^2 + 5^^^ - ^/ = 0. 



When the circumference of the source comes in contact with the caustics, the solu- 
tion of Eq. ([68]) has only one solution. This condition is 



where Or satisfies 



9{0c) = 0, 

dm 

de 



(69) 
(70) 
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X 
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1 
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0.2 


1.2 

1 
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0.4 
0.2 


0.2 0.4 0.6 0.2 0.4 0.6 
77 77 

Fig. 11. Boundary of the type I, III and V on the r] — 5j/io plane. The dashed curves satisfy 
= 0, the dot-dashed curves satisfy = 0, which are obtained by solving Eq. I I62I I of the 
sign + and — , respectively, while the solid curves plot Eq. I I65I I. In order to show the validity of 
the approximate approach, we plot the points of the critical configuration type II, the crosses are 
the exact approach (a), the triangles are the perturbative approach (b), and the squares are the 
approximate approach(c), respectively. The five asterisks in the upper left panel correspond to the 
value of r\ and 5j/io of Figs. 1 ~ 5 from the right to left, respectively. 



The solution of Eq. ^ is 



7c± 



COS 



^ f(52/io±\/32772 + 5yio^ 



then, Ea. (|5^ gives 



128772 



(71) 



(72) 
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This is the condition that the circumference of the source comes in contact with 
the caustics. Note that Oc± = dmin±, and g{Oc±) = Aj_. This means that the 
critical configuration type II (plus sign) and type IV (minus sign) appear when the 
circumference of the source comes in contact with the caustics. 

This condition can be transformed to the following relation between the critical 
line and the image. From Eq. (|46p . the central line of the image can be defined by 



6x 



(1 - a|0o(^)) 



dsS^{x,9) + Sy 10 cos ^ . (73) 



The critical line is defined by Eq. (|52p , then an intersection point of the critical line 
and the central line of the image satisfies 

(5yio cos 61 = — 2- . (74) 

OW^ x=l 

Within the approximate approach (c), using Eq. (jC.4[) . this condition gives 

SyioCosO = 2-qcos2e, (75) 

which can be solved easily, 



cos 6lcri± = ^ \ 5yw ± V 32772 + 5ywj ■ (76) 

Note that 6'cri± = 6'min±(= Oc±). 

The above behaviour of the critical configuration is obtained using the approxi- 
mate approach (c), but holds in the exact approach in a similar way. Indeed, these 
critical behaviour can be seen in Figs. [2] and H) These facts also guarantee the useful- 
ness of the approximate approach to investigate the lensing phenomena in a simple 
analytic way. 



4.3. Application of approximate approach 

In this subsection, let us summarise a few useful consequences, which are obtained 
using the approximate approach in an analytic manner. 

First, we consider the width of lensed images. Within the perturbative approach, 
the width of the image is 



From Eq. (j59p . the maximum width at 9 ^ ^max is 



(77) 



(1 - 5|0o(a;)) 

From Eq. ()61|) . the same maximum width appears at 6 = and tt. 
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Second, let us consider the angular size of the arc, which can be obtained by 
solving A^(6') — under the condition > 0, because A^(0) > is necessary for 
the appearance of the image. Using Eq. ([55]) . A'^{x — 1,9) =0 reduces to 

(5^/ - (?7sin26'- %rosin6l)2 = 0, (79) 

which can be solved easily with a suitable method. 

Third, we focus on the magnification o f th e lensed image. The magnification 
factor of an extended source is written (e.g.,^!^, 

_ J d^yl{y)fi{y) 

where I{y) is the surface brightness at y. In the case the surface brightness is a 
constant, i.e., — Iq, we may write 

Me = ■ (.^^i 

From the definition of the magnification /i, we have 

where we defined 

— = / _ r'dr' = -^(fc+ - 5x^). (83) 

Within the perturbative approach, we have 

.'^(^ fl\ 

(84) 



dA _ 1 2^A2(J,( 



dO ^^^^{l-dlUx)) 

In the approximate method (c) , we may substitute the expression ([55)1 into Eq. 

Fig. [T^ shows dA/dO as a function of 0. In each panel, a different set of the 

parameters Srg, Sy iq, rj and uq is adopted, as shown therein. The solid curves show 

the exact approach (a), the dashed curves show the perturbative approach (b), 

and the long dash-dotted curves show the approximate approach (c). The top left 

panel adopts the same parameters as those of Fig. [T] In the top right panel, 5rs is 

increased compared with the top left panel. This increase of Svg changes the four 

lensed images to one arc and the other separated image. In the middle left (right) 

panel, compared with the top left panel, uq is decreased (increased), by which the 

amplification is increased (decreased). In the lower left panel, compared with the 

top left panel, Syi^ is increased, by which the four separate images change to one 

arc and the other separated image. In the lower right panel, 77 is increased. 
. 2 

From Eqs. (f59| and (|6T1) . we have /S.^(x = 1,9) = Srs for 9 = 0, 0,„ax, tt, where 
the width of the image becomes maximum. Then, the maximum value of dA/d9 is 

dA 2 , , 
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Thus, dA/dO takes the same value at = 0, ^max, tt, which can be also seen from 
Fig.Hl 



20 
16 

dA 12 

8 



O 




20 |- 5r^ = 0.07, 5yio = 0.0 
77 = 0.15, Uo = 0.3 




T 



5r^ = 0.07, (5yjo = 0.09- 
77 = 0.15, Uo = 0.5 




- <5r3 = 0.07, t5yi„ = 0.1 
77 = 0.15, Uo = 0.5 




5r^=0.12,' (5yio=0 09- 
77 = 0.15, Uo = 0.5 



c5r,= 0.07, 6y,^ 
77 = 0.15, Uo = 0. 



= 0.09H 
9 



6r,= O.OV, <5y,„ 
77 = 0.25, Uo = 0- 



= 0.09H 
5 



7T 

e 



^7T O JL. 

2 2 



2 



e 



Fig. 12. dA/dB as a function of 6. A different set of the parameters Srs, Syio, V ^iid '"0 
adopted for each panel as shown therein. The solid curves are the exact approach (a), the dashed 
curves are the perturbative approach (b), and the long dash-dotted curves arc the approximate 
approach (c). 



4.4. Point source limit 

Finally, in this section, we consider the limit of a point source, which is given by 
imposing Srg = 0. In this limit, Eq. (|47p yields 



I ~ _ - — - — \ 2 
— dg6(j){x, 9) — Syio sinfl + 61/20 cos ( 

X 



(86) 
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For the existence of the solution, we must have A{x,9) = 0. In the approximate 
approach (c), with the use of Eq. (jC.3p . A{x = 1,9) =0 yields 



77 sin 29 — Syio sin 9 + Sy2o cos 9 — 0. 
This equation is equivalent to 



(87) 



... 2 — 2 

?7 sin 261 - V (52/10 + Sy2o sm{9 - 9o) ^ 0, 

where 9o = arctan((5y2o/<5yio) is regarded as an inclination angle when choosing 
the coordinate so that the point source is located on the j/i-axis. It is easy to 
solve Eq. ([55]) in a numerical method. When a solution of Eq. ([55]) is found, which 
we denote by 9p, from Eqs. (|46l) and ([50)1 . we obtain the image position and the 
magnification factor. 



Sx 



1 



d^5(j){x, 9) + Syio cos ( 



and 



(1 - 9|0o(^)) (5yiocos6' - 2?; cos 26* 



(89) 



(90) 



x=i,e=e 



respectively. 

For simplicity, we here consider the case ^j/20 = Oj that is, the inclination angle 
is zero. In this case, Eqs. ([5^1 and (lUHl) yield simple analytic expressions, as follows. 
Eq. dST]) reduces to 

?7sin26'-??/iosin6' = 0, (91) 
which gives the following solution to represent the angular position of the image, 



and 



2r; 



(92) 



(93) 



Eq. ([M)) gives the radial position of the image. For the solution of Ea. ([M| . we have 

' 2 



5x 



1 



dsScpix, 



Sy 



10 



2t] 



By substituting Eq. (l94l) into Eq. (|90| . we have the magnification factor. 



(l-9|0o(x)) 



277 



(94) 



(95) 



x=ifi=e 
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In the approximate approach (c), using Eqs. (|C.4[) and this magnification 

factor is expressed 



277(1 - 5|(/.o(5?)) 1 - 



6y 



10 



277 



(96) 



In the same way, for the solution of Eq. (1^51) , we obtain the radial position of the 
image, 

1 



'^^ (i-a|0o(s?)) 

and the magnification factor, 



d^5(j){x, 0) ± Syio 



2r,{l-dlM^)) -1± 



277 



(97) 



(98) 



where dl<j)o{x)\^_^ is given by Eq. (|B.1[) . where the sign + and — correspond to the 
solution 0p = and 0p = tt, respectively, from Eq. 



5. Summary and Conclusions 

We studied the perturbative approach to the strong lensing system, which was 
extended by Alard, in both the analytic and numerical manners. We investigated the 
validity of the perturbative approach by comparing with the exact approach on the 
basis of the numerical method, focusing on the shape of the image, the magnification, 
the caustics, and the critical line. The perturbative approach works well in the 
case when the ellipticity of the lens potential 77 is small and the configuration of 
the source is close to the that of an Einstein ring. At a quantitative level, the 
perturbative approach is valid at the 10 percent level for SyiQ ^ 0.2 and 77 ^ 0.3. 
We also demonstrated that the lowest-order expansion in terms of 77 also works well, 
which enables us to investigate the lensing system in an analytic way. 

We investigated the critical behaviour of the lensed images, by demonstrating 
the phase diagram of the different configurations of four separated images (type I) , 
an arc and one separated image (type III), and one connected ring image (type V). 
The critical configuration of type II appears during the transition from the type 
I to the type III, while the type IV appears during the transition from the type 
III to the type V. We investigated how the critical behaviour depends on the lens 
ellipticity, the source position and the source radius. We also demonstrated how the 
appearance of the critical configuration III and V is related to the condition between 
the source configuration and the caustics. The condition of the critical configuration 
was investigated in an analytic manner using the lowest-order expansion of the 
ellipticity 77 of the elliptical lens potential. 

The perturbative approach with the lowest-order expansion with respect to 77 
is useful to find the simple formulas which characterises the lensing system in an 
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analytic manner. We derived the analytical formulas of the arc width and the magni- 
fication factor. In the point source limit, the simple formulas for the image position 
and magnification factor were obtained. These formulas can be easily solved, which 
gave the simple analytic expressions in the absence from the inclination angle. These 
results will be useful to understand the gravitational lensing phenomena. 

In a realistic situation in reconstructing a gravitational lensing system, sub- 
structures in the lens might have to be taken into account. In the reference,!!^ 
Alard considered how a substructure affects a lensed image in the perturbative ap- 
proach. Even a substructure with small mass could make a change in the caustics 
and the lensed image drastically. It is an interesting problem how one can deter- 
mine the gravitational lens potential including substructures simultaneously. Here, 

17 

there is potentially a lot of room for improvement. This issue is outside the scope 
of the present paper, but need to be elaborated for a precise reconstruction of a 
gravitational lens system. 
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Appendix A. Formulas in the exact approach 

In this appendix, we summarize formulas of the lens equation in the exact approach 
(a) without any approximation. We consider the lensed image of the circumference 
of the circular source, whose center is located at ((5?/io, <5?/2o)- The source's radius is 
6rs ■ The circumference of the circular source is parameterized as 



with the parameter ip in the range < (p < 27r. As the image is parameterized as 



yi = Syio + 6rs cos ip, 
m = Sy2o + Srs cos (p 



(A.l) 
(A.2) 



X2 = X sm 



xi = X cos 0, 



9. 



(A.3) 
(A.4) 



the lens equation is 



Syio + cos ip = X cos 9 — cos 9—(j)(x, 9) + 

ox 

Q ^ ^ 

Sy20 + Si~s sm(p — xsm9 — sin6'— (/)(a;, 9) — ■ 




X o9 



(A.5) 



(A.6) 
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which yield 



X = 5yio cos 9 + 5y2o sin 6 + TF;</'(a;, 9) 

ox 



and 



±\ Sr 



tan(0 — ip) 



id_ 

xde 

1 d 



(A.7) 



■d9 



x,0) ~ Syio sin + Sy2o cos 9 



d 



(A., 



5yiocos9 - Sy2o sin( 
For the elliptical NFW lens model adopted in this paper, the potential is written, 

(A.9) 



4k., 1 



Uq X 



= — 2" ~ log 2arctanh^ 



1 + S 



where we defined ^ = UQXy/l — rj cos 26. Here, the expression of the case S < 1 is 
presented. Eq. (IA.9I) gives 



d(j) _ 1_ ( , - 
dS:~ u^n x[°^ 2 ' VI -S2 



90 J] sin 20 

9^ ^ 



Uq (1 — r]cos29) 



log- 



arctanh 
2 



l + S 
arctanh 



1 + S 



9^ 



4/^ g 2 ~ . _ „ I 



- 1 



1 



(1-S2)3/: 



-arctanh 



^1. 

,,2 ~2 



1 + S 

"2 



loa 



l-S^ 



arctanh 



2(-l + 2S2) 
VI -S2 
/ 2r7(cos26' - t;) 



1 + S 



V(l-'?cos2e')2 ^2 (l-?7Cos26')(l-S2) 



ulx'^rf sin^ 20 

""2>, 



4?7(1 - S2)(cos 29-T]) + 2r;2s2 sin^ 26* 
(1 -77 cos 261)2 (l-S2)3/2 

The case S < 1 is given by the analytic continuation. 



arctanh 



1 + S 



(A.IO) 
(A.ll) 

(A.12) 

(A.13) 

(A.14) 
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Appendix B. Formulas in the perturbative approach 

We here summarise useful formulas in the perturbative approach (b). For the ellip- 
tical NFW lens potential, we have 



— , 



—i 

dx 



9^. 



4k, 



log 



Mo 



1 



u, 







, -1 + 2^^ 



arctanh 



1 - uo 



x=i ui 



Aks rj sin 29 



arctanh 



arctanh 



log- 



z arctanh 



1- 












2 




1- 


MoT 




MoT 



1 ~ mq 

1 + Mo 



2772u§T2sin26' 



arctanh 



2"f 2 



4(1 - - r;2) 
ri'ulT^ sin^ 20 



1 ~ mqT 
1 + MoT 



m2T2 



arctanh 



MoT \ 
MoT/ 



1 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



where we defined T = ^1 — rj cos 29. 



Appendix C. Formulas in the approximate approach 

In the approximate approach (c), the potential is approximated at the lowest order 
of expansion with respect to rj. The following formulas are useful. By expanding the 
potential 5(f), Eq. (|45l) . we have 



5(t)[x,9) ~ ^77Cos2e'(-ilog^ - ^ arctanh 

M5 V 2 2 (1-m2)i/2 




(C.l) 
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at the lowest order of ry. Similarly, Eqs. (IB.2I) ^ (IB.4[) lead to 



dx 



rism26, 
- 2?7Cos26l. 



2(1-^2) {l-ulff- 



-arctanh 



1 - Wo 



l + uoj' 



(C.2) 
(C.3) 

(C.4) 
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